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Abstract 

We study the time evolution of different Laponite suspensions from a low-viscosity ergodic 
state to a viscoelastic non-ergodic state over a wide range of volume fractions and salt contents. 
We find that the evolution of non-ergodicity parameter (Debye- Waller factor) splits into two 
branches for all the samples, which correspond to two distinct dynamically arrested states. At 
moderately high salt concentrations, on the other hand, a third and new nonergodic state appears 
that is different from the above two nonergodic states. Measurement of the conductivity of 
i-rt ■ Laponite solutions in pure water shows that the contribution of counterions in the ionic strength 

is considerable and their role should be taken into account in interpretations of aging dynamics 
and the phase diagram. Based on these data and available data in the literature, we propose a 
new (non-equilibrium) phase diagram for Laponite suspensions. 
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I. INTRODUCTION 

Understanding the phase behavior of clay suspensions is of important technical and sci- 
entific interest. Characterization of different ordered and disordered phases formed by clays 
is of direct importance for various industrial applications such as soil mechanics and for 
instance the control of viscoelastic properties of materials with clay additives. Of more fun- 
damental importance is the study of the underlying mechanism behind gelation and glass 
formation that are both observed in clays. This potentially provides us with a deeper un- 
derstanding of dynamically arrested states of matter. Clay suspensions can typically be 
modelled as charged anisotropic particles such as disks immersed in an electrolyte, thus in- 
teracting via excluded volume, long-range electrostatic repulsions and weak (van der Waals) 
attractions. The phase diagram of anisotropic charged colloids such as clays and understand- 
ing the aggregation, gelation and glass formation appearing in such systems is a matter of 
considerable debate |l|, |2j, y, U, la, la, L3, la] • The specific clay system we study here, Laponite 
(a synthetic clay [9|) has been the subject of intensive study over the past decade or so. 

Laponite consists of crystalline platelets with a thickness of 1 nm and an average diameter 
about 30 nm and a bulk density of 2.6 g/cm 3 J9j. Each Laponite particle is a three-layer 
silicate composed of a central octahedrally coordinated magnesium-oxygen-hydroxide sand- 
wiched between two tetrahedrally coordinated silica-oxygen sheets. Isomorphic substitutions 
of the divalent magnesium atoms in the central layer by monovalent Lithium atoms lead to 
the formation of negative charges within the lattice, which is balanced by the sodium ions 
located at the surface. When Laponite is dispersed in water or any polar liquid, the polar 
molecules penetrate between interleaf regions, dissolving the interleaf cations and separat- 
ing platelet surfaces by hydration and electrostatic forces. Thus in the final suspension the 
Laponite surface has negative charge on the order a few thousand electron charge (in wa- 
ter), while its edges (depending on the pH) may have a small localized positive or negative 
charge generated by desorption or absorption of hydroxyl group where the crystal structure 
terminates. 

For a range of Laponite concentrations and salt contents, the dispersion of Laponite in 
water is followed by spontaneous evolution from an liquid-like state to a non-ergodic solid- 
like state. This process is called aging, meaning that the physical observables of the system 
such as diffusion of the particles and dispersion viscosity evolve with time. 

The aging dynamics during the fluid-solid transition in Laponite sus pen sions has been 
independently studied by manvgroups now @, |, |, 0, M, EH, El, El, H El, El, E3, El, El, 



20l . |2lL |22| . |23| . |2J, |25j, |26j, |27|, |28( using light scattering and rheology techniques. Perhaps 
the earliest studies on the phase behavior come from Mourchid et al. [2j and Kroon et al. 
ll| . Mourchid et al. attempted to characterize the phase diagram of Laponite suspensions 
based on rheological measurements [2|. They varied both particle concentration and ionic 
strength and by performing oscillatory shear measurements on samples one week after their 
preparation they defined a sol-gel transition line where the zero frequency elastic shear 
modulus increases remarkably. Their phase diagram gives a general overview, but we have 
come to realize that their method is flawed as the measurements were done after some 
arbitrary waiting time t w and viscoelastic properties depend on t w . 

On the other hand, Kroon et al. [ll( studied the aging of Laponite using dynamic light 
scattering (DLS) experiments. They measured a range of sample concentrations between 2.2 
and 3.5 wt% and found that all the samples show a similar aging behavior and evolve from 
an initially ergodic state to a non-ergodic state around a certain time (ergodicity-breaking 



point detected by changes in the moments of scattered intensity distribution) that decreases 
exponentially with increasing concentration. For a 3 wt% sample of Laponite, they reported 
the growth of the non-ergodicity parameter (fraction of frozen-in density fluctuations) from 
almost zero to approximately 0.8 at the late stage of aging. Bonn et al. |3j suggested 
that aging in samples with no added salt is due to strong electrostatic repulsions, leading 
to the formation of a low volume fraction Wigner glass. They determined the liquid-glass 
transition volume fraction as a function of ionic strength, assuming that the effective volume 
per particle can be estimated as ttR 2 Id (here R ~ 15nm is the particle radius and Id is the 
Debye length)by considering the Debye length as the particle thickness. If this exceeds the 
volume available per particle, which is estimated as irR 2 h/(f),(h = 1 nm the actual thickness 
<f) is the volume fraction) no free volume is available and thus a glassy state may emerge. 
The volume fraction for the liquid-glass transition is given by <p e ff ~ 0.5 and thus varies as 
dig ~ h/lf) oc J 1 / 2 where / is the ionic strength. Following the suggestion of Bonn et al. 
15 1, Levitz et al. 13j, by deionizing Laponite indeed found the evidence for a Wigner glass 



at very low ionic strengths consistent with the proposal of Bonn et al.. 

Nicolai and Cocard |l8|, ll9( have studied the aging at low concentrations of Laponite 



with added salt in the ergodic regime of aging. They observed that the scattered intensity 
increased with the aging time. At the late stages of aging, when the scattered intensity has 
become constant, the structure factor S(q) shows a power-law q-dependence characteristic 
of fractal structures for S(q). In a later paper the same group proposed a revised state 
diagram for Laponite suspensions based on visual observations and waiting time- dependent 
static light scattering experiments [4]. According to this study, the transition from liquid 
to "solid" (defined as the state that does not flow when tube containing the sample is 
tilted) occurs at much lower concentrations than what was proposed in the phase diagram 
of Mourchid et al. [2]. From their observations they argue that the origin of aging for all 
their measured samples with salt and without salt {C < 2 wt%) is gelation rather than glass 
formation [4|. 

The systematic study of Ruzicka et al. 22j, |23j was the first one to show that non-ergodic 



states of Laponite can exist at very low concentrations (C ~ 0.3 wt%) of Laponite, in 
the region which was proposed to be a sol according to the phase diagram of Mourchid 
et al. |2J. Using DLS, Ruzicka et al. 22j, |23j systematically studied the aging dynamics 



of both low and high concentrations and also varied the salt concentration. They showed 
that the intensity correlation functions at low and high concentrations evolve in a distinctly 
different manner and two distinct master curves have been identified. They suggested that 
aging at low concentrations proceeds by formation of a Wigner glass made of Laponite 
clusters, while at higher concentrations a glass is formed whose basic unit is a single Laponite 
particle [22] . They show that this result is not affected by the presence of salt and one still 



finds two distinct routes of aging 23]. However, based on their recent small angle X-ray 
scattering measurements, Ruzicka et al. [8j conclude rather that the arrested state at low 
concentrations should be called a gel and at high concentrations an 'attractive glass'. 

The study of Ruzicka et al. is the most comprehensive and systematic one up to now. 
However, their experiments considered only the ergodic regime of aging and they do not 
present any results in the interesting range of concentrations 1.5 < C < 2.2, precisely the 
range between 'dilute' and 'concentrated' systems. In a recent paper [7|, for samples in 
the range 0.1 < C < 3.6 we showed that in the non-ergodic regime the non-ergodicity 
parameter (Debye- Waller factor) also falls onto either of two distinct master curves. This is 
consistent with the classification obtained based on evolution of dynamic structure factor in 



the ergodic regime of aging [22j.The evolution of non-ergodicity parameter provides us with 
some valuable information about the nature of non-ergodic states. Using this information 
altogether with other data such as short-time diffusion and structure factor at low q limit, 
we identified the two observed distinct states as gels and glasses [7J. Furthermore we showed 
that for a range of intermediate colloid concentrations 1.1 < C < 2.4, the transition to non- 
ergodicity can occur in either direction (gel or glass). The distinction between glass and gel 
was mainly made on the basis of (i) the difference in the dynamic structure factor, notably 
the absence (gel) or presence (glass) of cage-rattling motion; (ii) the difference in static 
structure factor: S(q) showed power-law behavior as a function of q and a clear increase in 
time for the gel state. Both are consistent with formation of some sort of structure (cluster 
or network) in the gel, and are absent in the glassy state. 

It may be evident from these observations that the nature of non-equilibrium phases 
formed by Laponite suspensions remains ambiguous. Both gelation [2j, |4| and glass formation 
[3j, lla, [22j, [28| have been proposed to account for the aging process. Gelation and the glass 



transition have important similarities. Both are ergodic to non-ergodic transitions that 
are kinetic, rather than thermodynamic in origin, and distinguishing between these two 
types of non-ergodic states experimentally is a longstanding controversy [2|, y, U, la, la]. 
Here we will show that at least part of the confusion about glassy or gel-like behavior of 
Laponite suspensions finds its origin in that each group has only studied a specific range of 
concentrations or salt content. Furthermore, most of the studies have been performed in the 
ergodic regime of aging. Other serious discrepancies between the results of different groups 
arise from the fact that some of the measurements for determining phase diagram ignored 



the aging features of Laponite and were done after some arbitrary waiting time t w [2|, [13 
Another important confusion regarding the phase diagram in the literature is the neglect of 
the ionic strength that originates from the release of sodium ions from the platelets when 
Laponite is dissolved in water [3J. The ionic strength resulting from the release of sodium 
ions can be estimated from conductivity measurements. We will show in the following that 
the ionic strength of Laponite suspensions in pure water is relatively large and cannot be 
ignored for determining the 'effective particle size' from adding the Debye length to the 
size. Furthermore, different groups have used different grades (XLG and RD) and batches 
of Laponite which also can change the results slightly. Therefore this has lead to apparent 
contradictions for the results reported by different groups, that we attempt to clear up in 
the present paper. 

Here we have extended our previous study [7] to Laponite samples with added salt, 
since an important part of the discussion about the phase diagram is related to samples to 
which salt is added to screen the electrostatic repulsion between the Laponite particles. We 
report here extensive light scattering measurements during the evolution from an initially 
ergodic liquid-like state to a non-ergodic solid-like state on these samples. We show that in 
the presence of moderate amounts of added saltan addition to the two distinct nonergodic 
states (A and B) of Laponite, reported before [7|, [22f, even a third option (C) appears to 
exist for the system. 

We discuss the nature of the three distinct dynamically arrested states, A-C, and identify 
these, with the help also of data of other groups in the existing literature, in different 
regions of the phase diagram. This allows us to propose a unifying phase diagram for the 
non-equilibrium states of Laponite. 



II. EXPERIMENTAL 

We study charged colloidal disks of Laponite XLG, with an average radius of 15 nm and 
1 nm thickness. Laponite can absorb water, increasing its weight up to 20%. Therefore, we 
first dried it in an oven at 100°C for one week and subsequently stored it in a desiccator. 

We prepare a number of Laponite samples with different concentrations and salt contents. 
Laponite solutions without added salt are prepared in ultra pure Millipore water (18.2 MQ 
cm -1 ) and are stirred vigorously by a magnetic stirrer for one hour and a half to make sure 
that the Laponite particles are fully dispersed. The dispersions are filtered using Millipore 
Millex AA 0.8 /^m filter units to obtain a reproducible initial state [3j. This instant defines 
the zero of waiting time, t w = 0. 

The samples with added salt (NaCl, from Sigma) are prepared by diluting the Laponite 
suspensions in pure water with a more concentrated salt solution [l8[]. For instance, a sample 
of 0.8 wt %, 6 niM NaCl is prepared by mixing equal volumes of 1.6 wt% Laponite solution 
in pure water with the same volume of a 12mM salt solution. 

A standard dynamic light scattering setup (A = 632.8 nm) with a coherence factor close 
to 1 (^ 0.98) measures the time-averaged intensity correlation functions (Eq. ([1])) in VV 
mode, i.e. when polarization of incident light and scattered light are both perpendicular 
(vertical) relative to the scattering plane. 

< I(q,t)I(q,0) > t 

9 ^ t] = </(,,o)>? (1) 

where ()t stands for the time average. In the ergodic regime of aging this is related to the 
electric field correlation function i.e., intermediate scattering function f(q,t) through the 
Siegert relation gt(q,t) = 1 + P\f(q, t)\ 2 , where (3 is an experimental factor close to one 
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29] . For all the aging samples, there is a certain aging time after which the time-averaged 
correlation functions are no longer equal to their ensemble-averaged values, i.e. they change 
from one position to another in the sample. This defines the ergodicity-breaking point £ e &. 
This point is experimentally determined as the waiting time for which the time-averaged 
normalized second moment < I{t) 2 > t / < I(t) > 2 of the scattered intensity I(t) is not 
equal to one anymore, in other words the measured normalized correlation function g t does 
not decay from 1 to anymore [ll| . 

For aging times t a > t eb , we calculate the ensemble- averaged electric field correlation 
function i. e. intermediate scattering function f(q,t) from the time-averaged intensity 
correlation function gt{q,t) and ensemble-averaged intensity Ie measured by rotating the 
sample at different heights [30]. 

f(q,t) = 1 + {I t /I E ){[g t {q,t) - g t {q,0) + 1] 1/2 - 1}. (2) 

The measurements are performed at scattering wave vector q = ^p sin(|), where 9 = 90° 
is the scattering angle. The correlation functions are measured at a rate depending on the 
speed of aging of different Laponite suspensions. 

III. RESULTS 

Measuring the intensity correlations of scattered light from a large number of aging 
Laponite suspensions, one always observes two regimes of aging in the evolution of the 
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FIG. 1: Evolution of intermediate scattering function f(g,t) for three Laponite samples with salt 
at scattering angle of 90°. The symbols present the measured correlation functions at increasing 
waiting times (from left to right) that are (t w = 0.075, 5.7, 7.3, 8.8, 9.7, 11.9, 15, 19 and 500 days) 
for 0.8 wt%, 3 mM, (t w =11, 104, 153, 205, 255, 366 and 2854 min) for 2.5 wt%, 2 mM salt and 
(t w =9,44, 66, 90, 119, 164, 311 and 3900 min) for 0.8 wt%, 7.5 mM salt. In all panels, the lines, 
on the curves that decay to zero (ergodic stage), show the fits with f(q,t) = Aexp(— i/ri) + (1 — 
A)exp(-(i/r 2 ^). 



intensity correlation functions. In the first regime the system is ergodic, whereas the second 
regime corresponds to a non-ergodic (arrested) state. The cross-over from the former to the 
latter is visible in the experiments: the time-averaged normalized correlation function no 
longer varies between one and zero, i.e., a part of the degrees of freedom are frozen in on the 
time scale of the measurement. The waiting time for which the time-averaged correlation 
functions are no longer equal to their ensemble-averaged values (i.e., their values change 
from one position to another in the sample) defines the ergodicity-breaking time t e \,. Figure 
[Dshows the evolution of ensemble-averaged intermediate scattering functions f(q, t) for three 
different samples. In all the cases, the correlation functions evolve from an ergodic state to 
a non-ergodic state as the system ages. One can also observe that this generic behavior is 
not affected by the presence of salt. 

The intermediate scattering functions in the ergodic regime can be fitted with the func- 
tional form f(q, t) = Aexp(—t/ri) + (1 — A) exp(— (t/r?)^), in which T\ and r 2 represent the 
fast and slow relaxation times respectively [14j, |22J. In the non-ergodic regime, the aging 
rate of the system can be quantified by measuring the time evolution of the non-ergodicity 
parameter f(q, oo,t w ) = Hindoo f(q,t w ) [30|. In the absence of salt we have already shown 
that the evolution of non-ergodicity parameter in a range of samples with different Laponite 
concentrations collapses onto two distinct master curves when plotted as a function of re- 
duced waiting time (t w /t e b — l)(Fig. 2a). These branches were interpreted [7j as belonging 
to a colloidal gel and colloidal glass state based on different aging behavior in other mea- 
sured quantities. Here we refer to the branches for low and high concentrations as A and 
B, respectively. The observed differences between A and B in the absence of salt can be 
summarized as follows. 

• The static structure S(q) of B changes very little with waiting time, while that 
of A evolves dramatically. This is due to formation of network-like structure or 
aggregation. This difference manifests itself in the evolution of scattered intensity 
with waiting time, which grows for A but is nearly constant for B. |7J] 




FIG. 2: The evolution of the non-ergodicity parameter f(q,oo,t w ) versus reduced waiting time 
tw/teb ~ 1 f° r different Laponite samples a) without salt b) with salt. The colloid concentrations 
and salt contents are shown in the legend. The samples can be divided into two groups according 
to the evolution of non-ergodicity parameters. In part b the dashed and solid lines correspond to 
the aging process of group A and B, respectively, which are obtained from smoothed averaging 
over the data of Laponite in pure water. As can be seen, for most of the samples with salt (group 
C) the non-ergodicity parameter deviates from the glass line at long waiting times. The inset of 
panel b shows the difference in non-ergodicity parameter between the glass (line, averaged over a 
large number of experiments) and the attractive glass (data points). 



The short-time diffusion of particles in B decreases only slightly while it drops signif- 
icantly in A during the ergodic to non-ergodic transition [7J (See also Fig. \$&). 

The slow relaxation time of B grows exponentially with waiting time, while that of A 
grows faster than exponentially [71. |22||. 

The distribution of relaxation times is different between A and B: A has a broad 
distribution, whereas B has a double-peaked broad distribution of relaxation times. 

0. 

The difference between A and B is perhaps best reflected in their q-dependence of 



their structure factor at low q values and for late aging times. While the structure 
factor of B is flat at low q- values, S(q) of A is q-dependent, indicating that a structure 
has been formed [7|. 

The non-ergodicity parameter for A grows at a faster rate than for B. While the non- 
ergodicity parameter for A asymptotically reaches one, the non-ergodicity parameter 
for B reaches an approximate value of 0.85 for late aging times, indicating that there 
is still some freedom for the particles to move. This is suggestive of the "cage rattling" 
icture of glassy dynamics, also in agreement with the distributions of relaxation times 



m 



As explained in |7|, in view of the above measurements, it is tempting to identify A with 
a gel, and B with a glassy phase. This will be discussed in detail below; we now first consider 
the effect of added salt. 
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FIG. 3: (a) The ergodicity-breaking time t e f, as a function of the salt content for a few Laponite 
concentrations. i e (, decreases with adding salt. The dotted lines here are just guidelines for the 
eyes. (b)The scattered intensity relative to the toluene intensity I r as a function of dimensionless 
scattering vector qR for different amounts of added salt (NaCl) at a concentration of 0.8 wt%, 
as shown in the legends. These data are taken a long time after the samples have become fully 
non-ergodic. 



It is clear that adding salt (NaCl) to a given concentration of Laponite accelerates the 
aging. Figure Eh shows the ergodicity-breaking time for the three samples as a function of 
salt concentration. The effect is tremendous: by adding a few mM of salt, t e b can decrease 
by 4 orders of magnitude, with a roughly exponential dependence of the ergodicity breaking 
time on salt concentration. We can see the change of the slope in Fig. 3(a) for 0.8 wt% 
around 5 mM salt concentration, which reflects the crossover from A to C. 

Furthermore, looking at the structure factor at the late stages of aging (when the scatter- 
ing intensity has been stabilized, 50£ e 6 < t w < 100£ e &) for different salt contents, we find that 
with increasing salt the intensity increases and the wave- vector dependence of I(q) becomes 
more pronounced (Fig. [3b). The observed change of the structure factor with an increase 
in salt is plausible. The more salt we add, the more we suppress the repulsive interactions, 
therefore the attractive interactions play a more dominant role, leading to formation of 
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denser clusters and a more heterogenous structure. This result is consistent with the results 
of Refs. |4J, |20j in which the dependence of final structure on salt content has been studied 



for several Laponite concentrations (0.1, 1 and 1.5 wt%). 

Figure [2b shows that in spite of the accelerated aging, the evolution of non-ergodicity 
parameter f(q, oo), versus scaled waiting time tyJt e b still splits into two branches similar to 
those observed before for samples without salt [7|. 

However, looking carefully at Fig. [2b, it is evident that some of the samples with salt 
deviate from the 'glass' (B) master curve (obtained from the data without salt) for longer 
waiting times. These samples seem to evolve faster than the glass for long waiting times t w > 
3t e b and the non-ergodicity parameter reaches values higher than for the glass, indicating a 
blocking of rescaled particle motion. Measurements performed on these samples at very long 
waiting times show that the non-ergodicity parameter of these samples asymptotically reach 
the value 1. Hence, it turns out that in the presence of salt, the story is more complicated 
than the scenario sketched without salt [7J. If we look, for instance, at the scattered intensity 
as a function of time (Fig. H]), we find that for most of the samples with salt the scattered 
intensity increases and concomitantly their diffusion coefficient decreases with waiting time. 
Both increase of intensity and decrease of short-time diffusion are in principle indicative of 
the building up of structure, and thus suggest that a gel forms [7]: A=gel. Comparing, 
however, with the 'master curves' for the non-ergodicity parameter, we find that the high- 
salt concentration samples (Lap 0.8 wt% with 5 and 7.5 mM) behave more like B, whereas 
the low-salt concentration samples (Lap 0.8 wt% with 1 and 3 mM) are on the master curve 
of A. All of the 1.5 wt% Laponite samples should be glassy also, according to the non- 
ergodicity parameter criterion; however at least for the 1.5 wt%, 7.5 mM sample also a clear 
increase in intensity is observed. 

Plotting the slow relaxation time T2 normalized to its value at t w « as a function of 
scaled waiting time t w /t e b, we find that T2/V0 for all the samples, with or without salt, splits 
into two branches. We can see that for most of the samples belonging to the B branch of 
non-ergodicity parameter the slow relaxation time grows exponentially while for most of 
the samples belonging to the A branch T2 grows faster than exponentially, similar to what 
observed for samples without salt [71 122| . 

Comparing between different quantities, an inconsistency appears, which is always the 
same one. Looking at the non-ergodicity parameter, all the samples at moderate salt always 
behave like B. However, some of their other measured quantities such as scattered intensity, 
short-time diffusion and slow relaxation time, consistently behave as if the sample were 
A. This situation is indeed quite different from the one without salt as described in J7|. 
We see no 'hesitations' of the samples between two states in the sense that a sample that 
starts evolving in one direction ends up in the other one. Rather, all individually measured 
quantities consistently show an evolution in one direction. The data therefore suggest that 
although the sample has some definite characteristics of B, the other characteristics are those 
of A. 

Hence, to summarize, the addition of salt introduces new patterns in the aging behavior 
in the sense that there are samples which share some of the properties of A and some of the 
features of B. We call this new set of samples C. With this classification, with added salt 
the samples Laponite 0.8 wt% with 5 and 7.5 mM salt and Laponite 1.5 wt%, with 3, 5 and 
7.5 mM belong to the new group C. The samples 2.5 wt% with 1 and 2 mM salt behave in 
all respects identically to the group B samples: they are glassy. Also samples 0.8 wt%, 1 
and 3 mM behave exactly like as group A samples, and therefore are gels. 
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FIG. 4: (a) The evolution of short-time translational diffusion normalized to its initial value (t w ~ 
0) as a function of t w /t e b. The solid and dashed lines show the B and A line, respectively, obtained 
from smoothed averaging over the data of Laponite in pure water, (b) Scattering intensity at 
scattering angle 90° as a function reduced waiting time. So as to focus on the effect of aging, we 
have normalized the intensity to its value at the beginning of aging. 



IV. A PHASE DIAGRAM FOR NON-EQUILIBRIUM STATES OF LAPONITE 



As we discussed in detail above, different groups have studied the phase behavior of 
Laponite suspensions [2|, |ll|, [l9|, [22] , without a consensus emerging. Here, we would like to 
show that one can get a consistent picture putting all different pieces of information from 
different groups, despite the fact that each group has used different batches of Laponite and 
sometime different grades (Laponite XLG and RD). To demonstrate this point, we have 
plotted in Fig. Ek the ergodicity-breaking time obtained from our measurements [3], from 
Kroon et al.'s experiments [11] as well as t^ 3 denned by Ruzicka et al. 22|. Ruzicka et 
al. obtained t™ from fitting the waiting-time-dependent mean relaxation time r m with the 



general form T m (t w ) = T exp(B ™ ). Following the proposal of Ruzicka et. al. [22 



we 



fitted the mean relaxation time from our data to the above form in order to determine t^. 
Comparing t^ obtained from our fits with t^ obtained directly from the experiments, it 
turned out that for lower concentrations t^ m t^ and for higher concentrations t e f, ~ 0.6^^. 
can also be interpreted as a characteristic time for the transition from fluid- 
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FIG. 5: The evolution of slow relaxation time as a function of the scaled waiting time t w /t e b- The 
solid and dashed lines show the B and A line, respectively, obtained from smoothed averaging over 
the data of Laponite in pure water. 



like to solid-like state. Figure [6^, clearly shows that despite the difference in aging speed 
which is due to different batches of Laponite, the concentration dependence of ergodicity- 
breaking time found by the different groups is very similar. The two main reasons for the 
observed differences are most likely presence of additional water in experiments of Kroon 
et al. [ll| (The Laponite was not dried), and the difference in salt impurities between the 
very pure XLG (our experiments) and RD (Kroon et al. [H| and Ruzicka et al. [22]). A 
new phase diagram for non-equilibrium states of Laponite based on our characterization is 
shown in Fig. \5}p. Our data suggest the existence of three distinct nonergodic states A-C 
as demonstrated above. The Ruzicka et al. data also reveal the existence of two different 
arrested states (called IG1 and IG2 in their paper) for Laponite samples in pure water and 
at low salt content. Indeed the Phase IG1 of Ruzicka et al. corresponds to what we call A 
and their phase IG2 to B. The consistency between their and our data becomes even clearer 
if we plot the concentration at which the transition from A to B occurs as a function of 
added salt, as depicted in Fig. As can be seen there is a fair agreement for location of A 
to B transition line obtained from our data and the Ruzicka et al. data. 

Note that in the phase diagram of Fig. EJd the coordinate (y-axis) is the amount of added 
salt, while the interparticle interactions between particles are controlled by the total number 
of ions in the solution, i.e., the counterions released from the surface of the Laponite particles 
plus the ions from the added salt. To get an idea about the number density of counterions 
from Laponite, we have measured the conductivity of Laponite solutions. Figure [8^, shows 
the conductivity of Laponite solutions in pure water as a function of concentration measured 
at early stages of aging. We also measured the conductivity values for later stages of aging, 
before the samples become solid-like. We observed only very small changes, at most an 
increase in the conductivity of 5% as a function of aging time was found. The measured 
conductivity is mainly due to the Na + counterions released from surface of Laponite particles. 
The contribution of OH~ ions released from the edges of Laponite particles is relatively 
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and t™ oc t e {, extracted from Ruzicka et al. work 22|] as a function of concentration of Laponite 
in pure water, (b) Our proposed phase diagram based on light scattering data for non-equilibrium 
states of Laponite. 



small. Neglecting this contribution, the number density of Na + ions n can be obtained from 
^Na = &Nal ' ^Na^-, where fj,N a is the mobility of Na + ions (//jv a = 5.19 x 10~ 8 m 2 s _1 V -1 |3l| ) 
and e is the electron charge. The number density of Laponite particles can be estimated 
as rii = $ m / PhVh where <3> m is the mass fraction of Laponite particles, pi and vl are the 
density and volume of an individual Laponite particle, (pl = 2.6 g/cm 3 is used here and vl 
is calculated assuming a disk radius of 15 nm and thickness of lnm ). This allows us to plot a 
new phase diagram (see Fig. [SJa) with y axis being the total ionic strength / = 1/2 J^ riiZf, 
in which rij is the number density of ion species i and Zi is the charge of the ion species 
i. The abscissa is the Laponite concentration. As can be seen the ionic strength resulting 
from counterions is considerable and cannot be ignored. 
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Therefore the total ionic strength in Laponite solutions in pure water is much higher than 
the 0.1 mM salt below which the Wigner glass was predicted [3J and observed by Levitz et 
al. [2J. Taking into account the contribution of counterions resolves the confusion about 
the absence or presence of Wigner glass in Laponite suspensions in pure water. Levitz 
et al. prepared Laponite solution of extremely low ionic strengths by deionizing Laponite 
suspensions and immersing until the desired ionic strength. As a result they could observe 
a low volume fraction solid-like state. Measuring the structure factor with ultra-small angle 
x-ray scattering, they correctly identified this phase as a Wigner glass driven mainly by 
long-range electrostatic repulsions. 
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V. DISCUSSION 

To summarize, we have shown that Laponite suspensions can form different types of non- 
ergodic states (A-C) upon changing concentration and salt content. We have shown that the 
evolution of the non-ergodicity parameter (Debye- Waller factor) falls into distinct branches 
for all Laponite and salt concentrations. 

Now we come to the most important question of what is the nature of states A-C. For 
samples without salt, there are the two branches, A and B. The key difference is in the 
evolution of the static structure factor and translational diffusion coefficient with waiting 
time. The static structure factor and short-time translational diffusion of B are independent 
of waiting time, while the same quantities are strong functions of the waiting time for A. In 
addition, the slow relaxation time in B grows exponentially with waiting time, whereas it 
grows faster than exponentially in A. 

In group B samples (high Laponite concentrations, no or little added salt), the spatial 
structure is homogenous. This group show similar aging patterns in the late stage as seen 
in hard sphere glasses [30] and simulations of Laponite in glassy state 32J. For instance, 



even in the latest stages of aging particles maintain their freedom of rattling in the cages 
formed by their neighboring particle, as evidenced by waiting time independent short-time 
diffusion and a non-ergodicity parameter less than 1, which never exceeds 0.85 even at the 
latest stage of aging. 

We suggest to call this group a repulsive glass in the same sense as the glass formed at high 
concentrations of hard spheres. Note that although both attractive and repulsive interactions 
are present in all ranges of Laponite and salt concentrations, we believe that attractions do 
not play a dominant role at these relatively high concentrations, as is evidenced by the 
homogenous structure of these suspensions. This analogy becomes clearer if we plot the 
ratio of average interparticle distance between Laponite particles d to particle size D = 
2R = 30 nm versus concentration (See Fig. [9]). d is estimated as (irFPh/cf)) 1 / 3 , where is 
the volume fraction of Laponite particles. Figure [9] clearly demonstrates that d is very much 
comparable to particle diameter D. Adding the Debye length to the particle size makes this 
correspondence even better. Our interpretation of this glassy state at high concentrations 
is a jammed state which is appearing at much lower volume fractions compared to spheres 
due to anisotropic shape of disks and their large excluded volume effects. Interestingly Fig. 
correctly pinpoints the onset of deviations from glassy behavior (C ~ 2 wt%). 

Recently Ruzicka et al. [8] assigned group B as attractive glass. With this assignment, 
the transition from gel to glass can be explained simply by that the increase in the volume 
fraction of the particles leads to the decrease in the void size in gel and eventually the void 
size decreases nearly to the particle size (attractive glass) . Then, however, this scenario can 
explain neither the absence of the slowing down of the single-particle diffusion for group 
B nor the absence of the increase in the scattering intensity since attractive glass should 
accompany the finite-time bond formation between particles and slight increase in the scat- 
tering intensity. From this consideration, we suggest that the scenario that group B is a 
repulsive glass is more plausible. 

The assignments of groups A and C are much more subtle. In group A samples (low 
concentrations, no or little added salt), the aging behavior is distinctly different from that 
of group B. A dramatic decrease of short-time diffusion and remarkable increase of scat- 
tered intensity with waiting time is observed. The final structure in such samples is very 
heterogenous, suggesting the formation of a structure in the form of cluster or gel network. 
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These samples are also characterized by a non-ergodicity parameter which reaches the value 
1 roughly at a time equal to twice the ergodicity breaking time. All these features are 
consistent with the formation of a gel of Laponite particles, due to attractions between the 
particles [7J. The assignment of group A as gel has been made by Jabbari-Farouji et al. 
(7J and more recently also by Ruzicka et al. [8j As shown in Fig. 9, it is obvious that in 
this range the repulsive interactions are not enough to stabilize a repulsive glass. Neverthe- 
less, it is not so clear why a gel formed by attractive interactions appears in the region of 
weaker screening of repulsive interactions than the glass. The anisotropic shape of Laponite 
together with its rim of opposite charges (or neutral) may result in the aggregation in the 
concentration range. The formation of gel in such a ver y d ilute regime may be a consequence 
of the competing attractive and repulsive interactions [33|, |34|. Thus, we conclude that A is 
indeed a gel. 

There is a third group of samples, group C samples (moderate Laponite concentrations 
and high salt content), which also show a heterogenous structure but their aging behavior 
shares some of the features of group B (glass) and some of the group A (gel) . For example 
their scattered intensity increases with waiting time (a characteristic of a gel) while their 
non-ergodicity parameter evolution shows a similar behavior to that of group B (a glass). 
However unlike group B, the non-ergodicity parameter of the group C does not saturate at 
a lower value than 1, but keeps on increasing to reach the value 1 asymptotically at very 
large t w > 10t eb . 

Thus, samples in group C share some of the features of the glass, and some of the gel. 
This could be due to the fact that particles aggregate to form clusters, the diffusion of 
which becomes hindered progressively as the clusters grow. It has indeed been proposed 
that such a 'cluster glass' exists |l|, [Z 



35, 36, 37 



for which the size of clusters grows in 
time; this in turn makes that their diffusion significantly slows down. Together with the 
small amplitude of the motion, the diffusion mode inside a cage may become more and more 
difficult to observe. This may explain why the non-ergodicity parameter reaches 1 at very 
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late waiting times. In suspensions of a similar type of clay (monmorinite), Schurtenberger 
and his coworkers |l| found 'cluster fluids' in the corresponding dilute region of the phase 
diagram. 

Our experiments show that particles diffusion is suppressed in both group A and C at late 
times. However, we may need to take special care when interpreting the information of Fig. 
4. When we add salt, it is expected that even at t w = (in our definition) some clusters may 
have already formed. For example, the scattering intensity at 90° increases more steeply 
for group A than for group C, but the final scattering intensity is lower for group A than 
for C (see Fig. 3b). This means that at t w — the samples are more heterogeneous for 
group C than for A. This may be natural on noting that the attractive interaction inducing 
aggregation is stronger for group C than for A. On the other hand, C has characters of both 
glass and gel, where the slowing down of the dynamics is both due to aggregation (as in 
a gel) and steric hindrance (as in a glass). The coexistence of these two characters may 
be quite naturally explained by phase-separation- induced dynamic arrest: phase separation 
leads to a dense region of Laponite, whose composition is high enough for the formation of 



attractive glass. Such a scenario was proposed for suspensions of uncharged colloids [38 
The strong heterogeneity of the state C is also compatible with this scenario. So we assign 
the state C as attractive glass formed by phase separation (arrested phase separation due 
to the formation of attractive glass). 

Phase separation requires the presence of attractive interactions. At this moment the 
nature of the attractive interactions between Laponite particles is unclear. Possible sources 
are the van der Waals interactions and the attractions between the positive charge on the 
rim and the negative charge on the surface of Laponite particles. Indeed, recent experiments 
39( | have shown evidence for a short-range attractive potential in the effective interaction 
energy. 

At this point it is worthwhile comparing the aging features of our 'attractive glass' induced 
by phase separation with other attraction-driven glassy systems such as the attractive glass 
formed in hard spheres with added short-ranged attractions 40j, |4l|, |42|. In hard sphere 



systems with added attractions, attractive glasses are formed at moderately high volume 
fractions of particles and for strong enough attractions. This can be achieved in experiments 
by adding polymers that cause a depletion interaction [40|, |4l|]. Light scattering studies in 



these systems have revealed significant differences between attractive and repulsive glasses 



in both their static and dynamic properties. Pham et al. 40j showed that upon increasing 
the attraction strength, entering the attractive glass region (for a fixed volume fraction 
of colloids), the peak position of the structure factor shifts to a higher q- value and its 
height slightly decreases. This shows that average interparticle distance is decreased upon 
increasing the attractions and particles bond in clusters, implying that the average number 
of nearest neighbors should decrease (leading to the decrease in peak height), and 'holes' 
open up which render the structure more heterogenous on a length scale of a few particles. 
The increased heterogeneity is reflected in a slight increase of the structure factor at low 
q values. A similar trend is observed in our data upon increasing the salt concentration 
that screens the repulsions; this should be equivalent to increasing the attractions in the 
colloid-polymer system. 

The differences in dynamics of attractive and repulsive glasses is most evident in the short- 
time relaxations |4(J, |42J. The short-time dynamics of particles progressively departs from 



free diffusion upon increasing the attraction. In fact, for attractive glasses the particles are 
confined so tightly by attractive potential wells that short-time diffusion drops dramatically 
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compared to the repulsive glass at the same particle concentration 40]. This is consistent 
with our attractive glass samples for which a decrease of the short-time diffusion is observed 
(Fig. H]). Thus, our attractive glass shares some important similarities with attractive glasses 
in other systems. 



VI. CONCLUSION 

To summarize, we report that there are at least three distinct types of dynamically 
arrested states in Laponite suspensions. We specify the aging process towards these non- 
ergodic states in detail, using both static and dynamic light scattering. Our data indi- 
cate that the competition between short-range (van der Waals) attractions and long-range 
(electrostatic) repulsions in anisotropic Laponite particles leads to a rich non-ergodic state 
diagram and corresponding aging behavior. On the basis of our data on the static struc- 
ture factor and the dynamics of the aging, in conjunction with other observations in the 
literature, we propose that the three observed distinctly different arrested states should be 
identified as gel (A), repulsive glass (B) and attractive glass (C). 

The gel state is formed at low clay concentrations and low amounts of added salt. It has a 
spatially heterogenous structure as evidenced by our static light scattering measurements. In 
this case the aggregation of particles either in the form of a network-like structure or clusters 
is responsible for the aging process. The main characteristics of aging in a gel are dramatic 
slowing down of translational diffusion and a fast growth of non-ergodicity parameter to a 
fully non-ergodic state specified by non-ergodicity parameter of value 1. 

The glassy state forms in relatively high concentrations of Laponite and low amounts of 
added salt. Here, the aging dynamics of a glass has its origin in cage-diffusion process: for 
short times or small displacements normal Brownian motion is observed; however, for larger 
times or excursions, the particles are confined in effective cages formed by their neighbors. 
This becomes more and more difficult as time goes on, due to the fact that the system finds 
deeper and deeper energy minima during the aging process. On the other hand, even for late 
times the particles maintain their free rattling in the cage, as evidenced by a waiting-time 
independent short-time diffusion and non-ergodicity parameter different from 1 even at the 
latest stages of aging. 

Our study also suggests that a third non-ergodic state exists in Laponite suspensions, 
which we call attractive glass. This state is formed for moderately high amounts of salt 
shares some features of a glass and some of the gel. It has a heterogenous spatial structure 
similar to a gel while its dynamics is something between that of a gel and a glass. 

It is interesting at this point to discuss the relation of our light scattering measurements 
on these non-ergodic states with their rheological properties. Most interesting of course is 
to see whether a difference between the two 'glassy' states: attractive (B) and repulsive (C) 
glass can be found. As described in detail in Ref. [27J, we performed local microrheology 
(MR) experiments on samples belonging to groups B and C. In this technique, one looks 
at the (Brownian) motion of a probe particle, from which the visco-elastic properties of the 
surrounding medium can be inferred. It was found that that although the complex shear 
modulus shows a very similar frequency dependence for both types of samples, the local 
MR measurements reveal the differences in the structure. Local shear moduli obtained from 
samples of group B are independent of position in the sample while for a sample in group C, 
a significant heterogeneity in the sample develops as the sample ages. Therefore the shear 
moduli differ from one position to another in the sample. This provides one more piece of 
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evidence for the classification proposed here. 
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